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A notion of quantum multipole (in particular, dipole) noise is considered. Quantum 
dipole noise is an analogue of quantum white noise but it acts in a Fock space with 
indefinite metric. Quantum white noise describes the leading term in the stochastic limit 
approximation to quantum dynamics while quantum multipole noise describes corrections 
to the leading term. We obtain and study generalized quantum stochastic equations 
describing corrections to the stochastic limit which include quantum dipole noise. 
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1. Introduction 

There are numerous works devoted to the studying of the long time behavior in 
quantum theory. Recently an exact formula for the time dependence of certain 
matrix elements (ABC- formula) has been obtained in pQ. 

The leading term in the weak coupling - large time regime in this formula 
corresponds to the stochastic limit (see In the operator formalism the stochastic 
limit is described by the quantum white noise and corresponding quantum stochastic 
differential equations 0-0- 

To describe corrections to the stochastic limit one needs a new operator struc- 
ture. For free fields an appropriate new notion was introduced in [Jj. We call it 
the quantum dipole noise. In this paper we first study mathematical properties 
of the quantum dipole noise and then consider interacting case and show that the 
quantum dipole noise describes first order corrections to the stochastic limit. We 
introduce quantum multipole noise which describes corrections of any order to the 
stochastic limit. Moreover we introduce a new type of quantum stochastic differ- 
ential equations which involves quantum dipole noise. These quantum stochastic 
differential equations we call generalized quantum stochastic differential equations. 
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Let us remind that the quantum white noise is an operator valued distribu- 
tion with commutation relations proportional to Dirac's (^-function. The simplest 
example of the quantum white noise is given by the following commutation relations 

[b(t),b+(r)}=5(r-t). (1.1) 

The white noise operators are the subject of investigation in numerous works 

In the present article we study higher order corrections to the stochastic limit. 
We show that a new mathematical object: dipole noise operators - naturally arises 
as the second order terms in the expansion of collective operators a^ 1 (k,t/X 2 ) in the 
series with respect to the A. Here a ± (k,t) := e ±iw ^ k ' t a ± (k), where a~(k) and a + (k) 
are bosonic annihilation and creation operators. 

An operator valued distribution with commutation relations proportional to 
derivative of Dirac's 5-function has being introduced in [xj. We call it the quantum 
dipole noise. The simplest example of quantum dipole noise is defined by the 
commutation relations 

[c(^), C + (r)]=^ ( 5'(r-^). (1.2) 

We introduce and study also further generalizations of the white and dipole noise 
operators which we call multipole noise operators. They satisfy the commutation 
relations 

[c m (t),c+(r)}=6 nim eS^(r-t), (1.3) 

where m, n = 0, 1, co(t) = b(t) and c\{t) = c(t). 

White noise operators arise as the first order terms in the expansion of the free 
fields in the series with respect to A and given by the following limit (in the sense 
of convergence of correlators, see |2]) 

lim -a ± (k,t/\ 2 ) = b ± (k,t). (1.4) 

A — *0 A 

Here b ± (k, t) are the white noise operators. They satisfy the following commutation 
relations 

[b-(k,t),b + (p,T)] = 2TtS(uj(k))S(T-t)S{p-k). (1.5) 

Next terms of the expansion of the collective operators in the series with respect 
to the A have being found in [7j. They are the dipole noise operators c (fc, f) 
satisfying the following commutation relations 

[c~(k, t),c + {p, t)] = 2iri6'(uj(k))5'(T - t)S(p - k). (1.6) 

We will also use notations a(k) = a~(k), b(k,t) = b~(k,t) and c(k,t) = c~(k,t). 

In section 2 we will construct a representation of dipole noise H1.2|l in a pseudo- 
Hilbert space with indefinite metric. We will show that the quantum dipole noise 
is a well defined operator valued distribution in a Fock space with indefinite met- 
ric. Notice that spaces with indefinite metric are widely used in quantum field 
theory for quantization of electromagnetic and other gauge fields (Gupta-Bleuler 
formalism, [H], see also [HIHH]). 
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Using multipole (white and dipole) noise operators one can make the formal 
expansion 

a ± (k,t/\ 2 ) = Xb ± (k,t) + X 2 c ± (k,t) + o(X 2 ). (1.7) 
The base for this expansion is the following expansion (in the sense of distributions) 

1 e rf:r/A 2 _ 2ir5{x)8(t) + X 2 2mS'(x)5'(t) + o(A 2 ), 
A^ 

where x and t are real variables. 

Using expansion l|1.7|l of the collective (in this case free) operators we obtain 
the equation for the higher order corrections to the stochastic limit of the evolution 
operator U\{t) of system interacting with Bose field. Then we bring this equation 
to the normally ordered form which is convenient for calculation of corrections for 
matrix elements of the evolution operator. 

One has the expansion for the rescaled evolution operator into the series at the 
coupling constant A 

U x (t/X 2 ) = U (t) + XU 1 (t) + .... 

The normal form of the equation for the evolution operator Uo(t) in the stochas- 
tic limit is equivalent to the quantum stochastic differential equation in the sense 
of JT] and looks like 

dU (t) = (DdB+ - D + dB t - -fD+Ddt)U Q (t). (1.8) 

Here we use standard notations dB t and dB + for stochastic differentials of quantum 
Brownian motion; D and D + are bounded operators in system Hilbert space and 7 
is a complex number. 

We will obtain quantum differential equations of a new form which include quan- 
tum multipole noise. They will describe the higher order corrections to the stochas- 
tic limit. For example for the first order correction Ui(t) to the evolution operator 
one has the generalized quantum stochastic differential equation 

dUi(t) = (DdC+ - D+dC t )Uo(t) + (DdB+ - D + dB t - -fD + Ddt)U x {t). (1.9) 

We call it generalized quantum stochastic differential equation because it contains 
stochastic differentials of dipole noise which we denoted as dCt and dC+ . These 
operators are defined in a Fock space with indefinite metric. 

In Sect. 2 quantum dipole noise operators in the Fock space with indefinite metric 
are defined. An expansion of the free fields into the series containing multipole noise 
operators is obtained in Sect. 3. Dipole noise on the simplex is considered in Sect. 4 
wcorrections to the Green function are described in Sect. 5. Finally generalized 
quantum stochastic differential equations describing corrections to the stochastic 
limit are presented in Sect. 6-9. 

2. Quantum Dipole Noise in Fock Space with Indefinite Metric 

The commutation relations 1)1. 2f) for the dipole noise mean that one has a family of 
operators c g parameterized by complex valued functions of real argument g(t) from 
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some space of functions (for example from Schwartz space of functions £(K)) with 
commutation relations 

[c f ,c+]=<f,g>, (2.10) 

where 



<f,g>=ilf'(t)g(t)dt. (2.11) 

Notice that the inner product < ■, ■ > is non-positively defined, i.e. it is an indefinite 
inner product. In fact, for f(t) = u(t) + iv(t) with u(t) and v(t) are real valued 
functions one has 

</,/>= -2 / u'(t)v(t)dt. 



This quantity can be positive or negative. Therefore the natural space of represen- 
tation of algebra (|2.10() is a vector space with indefinite inner product. 

Let us describe a general construction of a Fock space with indefinite metric 
(see 0E31)- Let Tii be a Hilbert space over field C and (•, is the inner product 
in Tii ■ This inner product induces in Ti i a structure of complex Frechet space with 
Hilbert topology. An indefinite inner product on Tii is a hermitian scsquilincar 
form, i.e. a map Tii X Tii — > C with the properties 

< /, ag + 0h >= a < /, g > +0 < /, h > 



< f,9 >= <9,f> 

for all f,g,h 6 Tii and a, /3 € C. The form < •,• > is non necessarily positively 
defined. It is called an indefinite metric. We suppose that if F(-) is a continuous 
linear functional over the topological space Tii then it can be uniquely represented 
in the form 

F(g)=<f,g>, geTii (2.12) 

for some / € Tii. Moreover, we suppose that for any / £ Tii Eq. (|2.12() determines 
a linear continuous functional F(g) over Tii- If (|2.12|) holds then we say that the 
Hilbert topology is consistent with the indefinite inner product. In this case the 
complex Frechet space with Hilbert topology and with the indefinite metric is being 
called a pseudo-Hilbert space. The condition of 12.12fl is equivalent to existence 
of a bounded hermitian with respect to the inner product (•, ■)-h 1 linear operator 
r\ : Tii ^ Tii such that the inverse operator r/^ 1 exists, bounded and moreover rj 
satisfies the following property: 

< f,9 >= (/,1S)hh Vf,geTii. 

Notice that one can vary the inner product in the definition of Hilbert topology. 
In fact, two different inner products in Tii, sa Y ("j an d ("i ')'h > induce the same 
topology if 3 C, C > 0: 

C(f,f) ni < (f,f)' ni <C'{f,f)m- 

Moreover, one can choose the inner product in such a way that the operator rj 
satisfies an additional condition 

rj 2 = 1. (2.13) 
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We will denote in this paper as (v)wi the inner product for which operator rj 
satisfies (I2.13|) . 

Let H 1 s :— Sym Hi ® . . .®TL\ be a symmetric tensor product of n copies of Hi 
for n > 1 and C for n = 0. Operator of symmetrization acts on /„ = / W (g) . . . <g> 
as 

Sym f n = 1 V ® . . . ® =: /„ (2.14) 

7T 

with summation over all permutations of the set {1, ...,n}. The indefinite inner 
product of two vectors /„, g„ € 7^ 3 for n > is given by 

< fn,9n >■= {fn,r)® 9n) H ®S 

where rj® := rj ® ... ® r\ is the tensor product of n copies of rj. For n = 0we define 

< fo,9o >•= fo9o 
Then we can define the Boson Fock space with indefinite metric 

oo 
n=0 

We mean that any vector <f> G J- (Hi) is a sequence of the form 

4> = (/o,/i, •■■,/«,■■■) 

with the convergent series 

ll^ll 2 = x>»'^)<J <°°- 

n=0 

Here /o G C is a complex number, /„ given by (|2.14|) and is being called the n— 
particlc component of the vector (j). The indefinite inner product of two vectors 
<\> = {fo, A, ■ • ■) and -0 = (g , 5i, ■ ■ ■) is g iven by 

OO 

< (j>,iP >=^2 < f n ,g n > . 

n=0 

The series is convergent because < f n: gn >— {fn,V® 9n)n.! and ?7 2 = 1. 

Creation and annihilation operators and Cf act on n-particle component of 
the vector tj> by usual formulae 

(Cffn)n+i = Vn+l Sym / ® /„ 
1 - 

( C/ /„)n-i = -= V < /, /« > Sym /« ® . . . /« . . . ® / W. 
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Here means that is missed. One has the canonical commutation relations 
(CCR) 

[c/,c+] =</.5> 
and also on finite vectors one has the relation 



< Cf(f>,1p > — < <f>,Cftp 



> 



which means that the annihilator c/ is the adjoint of the creator Cj with respect to 
the inner product < • , • > . 

In order to construct a representation of algebra (|2.1l)[l let us choose Tti to be 
the completion of the Schwartz space S'(R) with respect to inner product 



(/>0)wi :=^f \r\f(r)g(r)dr. (2.15) 

Here 

f(r) = j e^f(x)dx 

is the Fourier transform of f(x). We introduce an indefinite inner product on Tli 
by the formula 

<f,9>=i JjW)9{t)dt (2.16) 
which can be rewritten in the form 

< f ,9> = h; J T ^ r )-9( T ) dr - 

In this case the operator rj acts on the Fourier transform of the function f(x) as 

ivf)(r) = sign(r)/(r) 

where sign(r) = 1 if r > and —1 in the opposite case. Clearly this operator 
satisfies H2.130 and all conditions which are necessary to guarantee (|2.12(l . There- 
fore the Hilbert topology which given by H2.15(l is consistent with the indefinite 
metric lj2~TH)l . 

Any vector (f> £ J- (Hi) is a sequence of (symmetric) functions, 

= (h, fl(tl), ■ • • , fn(h, ■ ■ -,t n ), ■ ■ .). 

Here /o € C is a complex number as above. 

According to general construction described above creation and annihilation 
dipole noise operators and Cf act on rt-particle component of the vector <f> by 
formulae 

^ n+l 

(Cf fn)n+l(h, ■ ■ ■ , t n +l) = , j = / J f(U)fn(t\i ■ ■ ■ , U, ■ ■ ■ , *n+l)- 

1 y/n + 1 r— ^ 

v 1 

(c// n )n-l(*l, ■ • - ytn-x) = i\/n J f'(t)f n (t,ti, . . . ,t n -i)dt 
Here U means that argument ti is missed. 
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Remark 2.1. Motivated by the investigation of the space with the sesquilincar 
form p.llfl let us note that it is interesting to consider the following approach to 
solution of the Dirichlet problem. Let us consider the Laplace equation 

- Au(x) = f(x), leficf 1 (2.17) 

where f2 is an arbitrary domain in R d . Denote H the Hilbert space obtained after the 
completion and taking the factor space of the space of test functions with compact 
support T>(Q) with the scalar product 



h 

We write Eq. (|2.17|) in the form 

= (/,¥>), f€V(Q) (2-18) 

where 



(f,<P) = J f{x)<p(x)dx. 

Proposition 2.1. Let f G Li oc (Q) and there exists a constant Co > such that 

|(,/»| 2 < C a (<p,(p)H V^D(!1). 
Then there exists a unique solution u € H. of Ea. $2.18\l . 



3. Multipole Noise 

In this section we will show that the quantum multipole noise naturally arises from 
studying of the higher order corrections to the stochastic limit. 

Let L 2 {R d ) be the Hilbert space (with inner product (•,•)) of square integrable 
on K d functions. Denote then 

oc 

the Boson Fock space and a(/), a + (g) are annihilation and creation operators in T 
with commutation relations 

[a(f),a+(g)] = (f,g). (3.19) 
Here /, g S L 2 (R d ) and we will use the notation 



a(f) = J f(k)a(k)dk, k e R d . 

For a given complex valued test function g{k) on M. d and for smooth real val- 
ued function uj(k) define the (smeared) time dependent creation and annihilation 
operators as 

A + {t) = [ g(fc)e^( fc ) t a + (fc)dfc, A{t) = [ g{k)e~ luj{k ^ a{k)dk. (3.20) 
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The test function g(k) in physical applications plays a role of formfactor describing 
structure of interaction. The function u>(k) has meaning of dispersion law so that 
these operators are free evolutions of smeared creation and annihilation operators. 
We assume that the surface ui(k) — is nondegenerate, i.e. gradw ^ on the 
surface uj(k) = 0. 

Let us consider the commutator of rescalled creation and annihilation operators 

A x (t) := \A(t/\i), A+(t):=jA+(t/X 2 ) 

(we call them collective operators), where A is a real parameter. One has 

[A x (t),A+(r)]=± J d%(fc)|V i£ ^. 

In order to investigate the behavior of this commutator for small A let us consider 
the sequence of functions on R 2 indexed by real number A ^ 

h(x t t) = ^e^/A 2 . 

For any A ^ f\(x, t) is locally integrablc. Therefore it defines a distribution from 
S"(R 2 ). In the sense of distributions one has the relation ([2]) 

limfx[x,t) = 2ir5(x)5(i). 
We will prove the following generalization of this relation. 

Theorem 3.1. One has the asymptotic expansion (in the sense of distributions) 

1 „«t/A 2 



J_ e «*/A 2 = 2n J2 ¥U-6M{x)8M(t). (3.21) 

n— 



The asymptotic series \3.21\l means that for any test functions f and 4> from the 
Schwartz space S(M) and for any N S N 



N '-\2\n 



1 ;e txt / x2 f(x)^(t)dxdt-2TT^^^f^(Q)^ n \0)^o(X 2N ' 



Proof. Let f(x) and 4>(t) be test functions from the Schwartz space S'(R). Then 
one has 

^e ixt / x2 f(x)<p(t)dxdt = J f(r)^(X 2 T)dT. 
Here 

f(r) = f e^f(x)dx 



is the Fourier transform of f(x). It is well-known (see that /(r) € S(M) and 
the inverse Fourier transform is given by the formula 

f(x) = ^ / e-^f(r)dr. 
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Since <fr(t) 6 S(M) one can write the Taylor expansion for </>(A 2 r) in the Lagrange 
form 

N <\ 2 -r\ n (\ 2 t) N+1 

*(A a r) = #0) + £ ^J-^(O) + ^^-^(^(r)) 

n— 1 

where £a(t) € (0, A 2 t). Using this expansion one gets 

f N \2n r 

/ /(r)0(A 2 r)dr - ]T — ^">(0) / r n f(r)dr 

- Ar+1 /(r)0( A ' +1 )(a(r))dr. (3.22) 



iV 

E 

n=0 

A 2(JV+1) 



(iV + 1)! ' 7 
Now notice that 

' ' T n f{T)dT = 2m n f (n) [0). 



Then LHS of l|3.22|l becomes equal to 



N (i\ 2 \ n 

re Mt/A f{x)(t>(t)dxdt - 2nJ2 ^y-/ (n) (0)<£ (n) (0) 



n=0 



It is clear that RHS of is o{X 2N ). □ 

Remark 3.2. Let us denote non smeared creation and annihilation collective op- 
erators as 

af(k,t) := ia ± (fc,t/A 2 ) 

so that 

A$(t)= [ g(k)a+(k,t)dk. 



The expansion (j3.21|l for f\(x, t), which can be symbolically rewritten in the form 



le^ =27rexp (a 2 ^)^), 
induces the corresponding expansion for the non smeared collective operators 

N 

a x (k,t)=J2\ n c n [k,t)+o(\ N ) (3.23) 

such that c n {k, t), n = 0, 1, . . . satisfy the commutation relations 

[c n (k, t),c+(p, r)] = 6 n , m ^-6^(u(k))S^(r - t)S(k - p). 

n\ 

Remark 3.3. It would be interesting to study in that sense one can speak about a 
spectral decomposition of initial Fock space T into a family of spaces corresponding 
to the multipole noise creation and annihilation operators. 
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We will use also the notations 

co(k, t) = b(k, t), c\(k,t) = c(k,t), C2(k,t) = d(k,t). 

Here b(k,t),b + (k,t) are the white noise operators: 

[b(k, t), b+{p, t)] = 2TrS{uj(k))5{r - t)5(k - p) 

and c(fc, t), c + (fc, t) are the dipole noise operators: 

[c(k,t),c + { P ,T)} = 2iri6'(uj(k))5'(T-t)5{k-p). 

The second order terms in the expansion l|3.23[) (we denote their as d(k,t) and 
d + (k,t)) satisfy the commutation relations 

[d(k,t),d + {p,r)] = -ir5"{u(k))6"(T-t)5(k-p), 

For the smeared time dependent operators one has an expansion 

N 

A(t/\ 2 ) = £ A" +1 c„(t) + o(X N+1 ). (3.24) 

n=0 

The operators c n (t) satisfy the commutation relations 

[c n (t),c+(T)}=6 1hm j n S ( - n \T-t) 

with complex numbers 

oo 

7 n = t^ J a n da I ' e iau ^\g(k)\ 2 dk. 

— oo 

The expansion (|3.24l) holds for any TV G N in the case then the set of station- 
ary points of function u)(k) has zero intersection with supp g(k), i.e in the case 
|g(fc)| 2 G S(fl) for any G M. d such that stationary points of ui(k) do not belong 
to J7. The set of stationary points of ui(k) is defined as JC :— {fco | gradw(fco) = 0}. 
The expansion holds for N = 1 and |.g(fc)| 2 G S{R d ) for d > 3. 

4. Dipole Noise on the Simplex 

Some modification of Theorem 13. II is useful for getting the differential equation for 
the first order correction. This modification is about distributions on the standard 
simplex (see 

Let us remind the construction for distributions on the standard simplex. Define 

Co := j> : K+ -> C | <p = a.e.} 

C\ := {(f> : K + — > C | 4> is bounded and left — continuous at any t > 0} 
C := linear span of {Co U Ci}. 
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For any a > define <5 + (- — a) as the unique linear extension of the map 

S + (- - a) : <f> G C\ — > 0(a) 

<*+(•- a) :</>eC o ^0. 

One has the following result: in the sense of distributions on the simplex there 
exists the limit (see for more discussions) 

a ^.-^ = _l_ j+(( _ o) . (4 . 25) 

In order to introduce the derivative 8' + {- — a) we will define C\ as 

n 

Ci := {(/) : M + -> C | 0(t) = V) 0i(t)0[ O , ai j(t) for some n G N, < a, < oo and 

i—l 

4>i(t) is bounded, with bounded and continuous first derivative ^(t)}. 
This definition means that for any G C\ and a > left-derivative 

<PlW •= llm n : 

exists. Moreover it is left-continuous and bounded. 

With such defined Cl we define C as linear span of {Co U Ci}. Clearly C C C. 
We will denote the restriction of 5+(- — a) to C again as <$+(■ — a). 

For any a > define <5^_(- — a) as the unique linear extension of the map 

S' + (- - a) : d> G C x - -^(a) 

£+(■ - a) : G C -► 0. 

With such defined <$+(• — a) the first order correction to l|4.25[l (proportional to 
A 2 ) is given by the 

Theorem 4.2. One has the following limit in the sense of distributions on the 
standard simplex (a > is a positive number) 

l im J_ ( i_ e «(*- Q )A 2 - — -S+(t - a)) = S'(t - a) — —^—. 

T/izs means that for arbitrary test functions f(x) G S(M) and <f>(t) G C one /ias the 
limit 



lim 

a^o A 2 



i | dx ( dtf{x)${ty x(t - a),X2 -{5+{--a),<f) J dyf{y) 



o 

-(#.(.- a), 0) / (4.26) 



-ffere 

/(y) = / e™yf(x)dx. 
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Proof. The Eq. (|4.26l) clearly holds for € Co- Therefore let us consider the 
case € C\. It is enough to consider of the form cj)(x) = (j>i(x)Q[ Qjai ](x) 
with (j>i bounded, with bounded and continuous first derivative. In this case 
(S + (- — a),0) = 0(a) and (S' + (- — a),0) = —cf)' L {a). Making change of the vari- 
able y = (t — a)/A 2 in the first integral one gets 



L := lim — 

A^O A 2 



dx 

oo -a/A 




2 , 



dyf(x)^y + a)e^y-cb(a) / dyf(y) 



lim 

a^o A 2 



o o 

dyf(y)<f>(\ 2 y + a) - 0(a) / dyf(y) 



-a/A 2 



Let us consider the case a > a\. In this case 0(a) — and function 0(A 2 ?/ + a) 
is equal to zero outside of the interval [—a/ A 2 , (a\ — a)/ A 2 ]. Therefore one has 



(a!-a)/A 2 



L = lim — 

A^O A 2 



dy/(y)^(A 2 y + a) = 



-a/A 2 



because of boundness of 0i and Lemma l4~Tl In this case <fi' L (a) — so RHS of (|4.26l) 
also equal to zero and Ij4.26(l holds. 

Now let us consider the case a < oi. In this case 0(a) = 0i(a) and one has 



L = lim — r 

A^O A 2 



ii 

dyf(y)M^y + a)~Ma) I dyf(y) 



-a/A 2 



lira —r 

A^O A 2 



-a/A 2 o 

*//(?/)<MA 2 ?/ + a) - / dyf(y)M^y + a)-Ma) I dyf(y) 



lim —r 

a^o A 2 



-a/A 2 

dyf(y) (0i (A 2 y + a) - 0! (a)) - f dyf{y)^ (X 2 y + a) 



Using Lemma f4. II we get 



L = lim — 

a^o A 2 



dyf(y) (0i(A 2 2/ + a)-0 1 (a)) 



Using Lagrange's theorem one gets 

1 (A 2 y + a)-0 1 (a) = A 2 #' 1 (£ A ( 2 /)) 
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with £\(y) € (a + A 2 y, a). Therefore 



u 

lim | dyyf{v)&{tZx(y)). 



Now since 0'i(?a(2/)) is bounded and /(y) G one can apply Lebesgue's domi- 

nated convergence theorem. Finally one has 

o o 
L = 4>[(a) [ yj{y)dy = <P' L {a) [ yj{y)dy.U 



Lemma 4.1. For arbitrary test function f(x) e S(R) and two positive numbers 
b > a > one has 

b/\ 2 



J dxf(x)=o(\ 2 ). 



a/A 2 

Proof. In fact for large enough x and Vn <G N 3C = Const such that |/(x)| < 
Let us choose n — 3. Then one has 



lim -J 

a^o A 2 



fc/A 2 



dxf(x) 



a I A 2 



6/A 2 6/A 2 

If C f dx 

< lim — r / cfal/fa;)! < lim — r / — r = 0. 

~ a^o A 2 7 ~ A^o A 2 7 a; 3 



a/ A 2 



a/A 2 



Remark 4.4. We get the following relation 



1 e «(t-a)/A 2 = tf (t _ a ) 1 + iA 2 £,(i - a) j- - 1 - + o(A 2 ) 
A 2 +V 'i(ar-iO) +V ; dxi(x-i0) V ; 

The corresponding expansion for the collective operators has the form 

A(t/\ 2 ) = Xb(t) + \ 2 c(t) + o(A 2 ) 
with causal commutation relations 

[b(t),b+(r)]= l0 S + (T-t) 

[c(t),c+(r)}= 7l 5' + (T-t). 
Here complex numbers j n are given by 

o 



In 



J o n da J e la ^\g{k)\ 2 dk. 



□ 
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5. Corrections to the Green's Functions 

The master field in the standard formulation of the stochastic limit is defined as 
the limit as A — > of the Wightman correlation functions: 

<0|Ai(tl)...^(<n)|0>. 

For some models this limit is trivial (equals to zero). However this does not mean 
that the stochastic limit of such models is trivial because we can consider the limit 
at A — * of the chronologically ordered correlation functions (Green's functions) 

<0\T(A\(t 1 )...Al(t n ))\0> . 

Here the T-product is defined as 

T{A\{t l )...Al{t n ))=A^{t ll )...A^{U n ) 

if tii > ■ ■ ■ > tin ■ 

For the non smeared collective operators aj~(k, t) from Theorem [O] one has the 
following expression for Wigthman functions 

<0\a x (p 7 t)a+(k,T)\0>=^e^- t ^/ x2 S(p-k) 

= 2tt5(t - t)6(Lu( P ))6{ P -k) + \ 2 2iri6'(T - t)6' '(w(p))6(p - k) + o(A 2 ). 
On the analogy with Theorem 13. II it can be proved that 

Therefore one has for Green's functions 

< 0|T (a x (p, t)a+(k, r)) |0 >= ^6(r - Qe^-W)/* S (p _ fc) 

- «(r - O^-^P - *) + AV(r - V^^SiP - *) + (A 2 )- 
The first term in this expansion was obtained in 

6. A System Interacting with a Reservoir 

In the last sections we derive a generalized quantum stochastic differential equation 
for corrections to the stochastic limit of the evolution operator U(t). We will apply 
the dipole noise operators described above for study of the higher order corrections 
to the stochastic limit. Then we will bring these equations to the normally or- 
dered form and rewrite it in the form of generalized quantum stochastic differential 
equations. 

The study of the evolution operator U(t) describing a system with interaction 
is very important in statistical physics and quantum field theory. 

One of the methods of investigation of behavior of quantum system interacting 
with a reservoir is the stochastic limit method (see for details). One considers 
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quantum models with small coupling constant A for large time t. The idea of the 
method is in systematic application of special limiting procedure t — > oo and A — * 
such that A 2 i = r = Const. The result of such limiting procedure is being called 
the X 2 t— limit or stochastic limit. The limiting equation for the evolution operator 
is a quantum stochastic differential equation. 

Using this limiting procedure one can study dynamics of real physical system 
with small but finite coupling constant for large time t ~ t/A 2 without perturbation 
theory. For some important models dynamics in the stochastic limit approximation 
became integrable. 

The first rigorous result about interaction of a system with a reservoir where 
the role of A 2 i-rescalling begun emerge is due to Bogoliubov |T2]- Friedrichs, in 
the context of the now well known Friedrichs model |14|. was lead to consider the 
scaling limit 

t — ► oo, A — ► 0, A 2 i = t = Const 

by second order perturbation theory. This rescaling was used to derive the master 
equation by Van Hove |T5|. Prigogine |T5] . 

The stochastic limit method leads not only to the master equation but also to 
the equation for the dynamics of the reservoir, see [21 1171 H%| . 

We will apply the dipole noise operators described above for study of the higher 
order corrections to the stochastic limit. We will consider the quantum model of 
iV-level system (atom) interacting with a boson field. 

Let us describe more concretely this model. Let T as above be the Boson Fock 
space over the one-particle Hilbert space L 2 (M. d ). Let C N be the Hilbert space of 
iV-dimensional complex vectors. Then 

$ = C N ® T 

is the Hilbert space of the model which we consider. 

The total Hamiltonian of the model H\ is (it acts in $): 

H x =H a + XV. (6.27) 

Here Hq = 1 £g) Hqh + Hqs <8> 1 is the free Hamiltonian with 

H 0R = J uj(k)a + (k)a(k)dk; k G R d (6.28) 

(d = 3 in the physical space), A is a coupling constant, and the interaction V is 
given by 

V = i{D®a + (g)- D+ ®a(g)). (6.29) 

Here the operators a + (g) and a(g) describe creation and annihilation of bosons; D 
and D + are N x N- matrices describing transitions between atom levels. The func- 
tion g(k) is a formfactor (complex valued test function) describing the interaction 
of the iV-level system with the reservoir. This Hamiltonian has being applied in 
quantum optics, tunnelling processes etc (see ^3 GDI for details). 

The function uj(k) is a dispersion law. We suppose that the equation u(k) = ujq 
determines regular surface in R d (oj e K will be defined below). For example 
ui(k) = Vfc 2 + m 2 for massive bosons (m is the mass of boson) etc. 
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Dynamics of the system with H\ is being determined by the evolution operator 
U x (t) = e ltH o e ~itH x _ It ig thc so i ut j on f the equation (see 0): 



U x (t) = l-iX J V(t')U x (t')dt'. (6.30) 
o 

Here V(t) = e ltH °Ve~ ltH ° is the interaction operator in the interaction picture. 

We will consider two particular types of the model described above. For the 
first type N is arbitrary natural number and the matrices D and D + are connected 
with H s by the following conditions for some ujo & R (the so called rotating wave 
approximation condition): 

e itHos £) e -itH 0S _ jj e -ituj 
g itH os JJ+ e ~- HHqs _ jj+ e itu>a 

Under these conditions one immediately has 

V{t) =i{D® A+{t) - D+ ® A{t)) (6.31) 
with the time dependent creation and annihilation operators 

A+(t) = { g{k)e l ^-^ t a+{k)dk, A(t) = / g(fc)e^W-^*a(/c)dfc. 



After time rescalling t — > t/X 2 these operators according to results of the sections 3 
and 4 have an asymptotic expansion 

A(t/\ 2 ) = Xb(t) + X 2 c(t) + o{\ 2 ) (6.32) 

with white and dipole noise satisfying the (causal) commutation relations which 
follow from Remark 14.41 These relations are 

[b(t),b+(r)}= lQ 6 + (T-t) (6.33) 

[c(t),c+(r)}= ll S' + (T-t) (6.34) 
and other commutators are equal to zero. Here 70 and 71 are complex numbers 



* = / dk J g(fc) ' 2 n = / d ° I dke^~^\g{k)\ 2 
J uj(k)-ujQ-iO J J 



7i 



dk r / J 9m2 , m2 =- / daa I dke^ k ^\g(kt 



— OO 

Usual commutation relations follow from Theorem l3.ll 

[b(t),b+(T)]=j 5(T-t) (6.35) 
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[c{t) lC + (T)]^%5'(r-t) (6.36) 

with 

oc 

70= f da f dke ta ^-^\g(k)\ 2 



71 = -/ daa \ dke la ^ k) -^\g{k)^ 



The second type of the model is the spin-boson model. This model corresponds 
to the case N = 2 and no rotating wave approximation condition is assumed. 
The spin-boson Hamiltonian in the simple but non-trivial case e = has the form 
(see [THE!]) 

H x = -\&Vx ®1 + 1®Hor + \<t z ® (a{g) + a+(g)). (6.37) 

Here the free Hamiltonian of the boson field Hqr is given by l|6.28|l , a x , a z are the 
Pauli's matrices, and A > is a positive number. This model is widely used in 
physics and chemistry and for example describes a dynamical model of two-level 
system coupled to an environment (see |17l I21j for example) . 

For both these models we study asymptotical behavior of the evolution operator 
U\(t/X 2 ) with rescalled time. So we study the expansion 

U(\,t) := U x (t/X 2 ) = Uo + XU^ + oiX 2 ). (6.38) 

It was established J2 that Uq satisfies the quantum stochastic differential equa- 
tion. White noise operators play an important role in this consideration. In the 
present article we show that higher order correction to Uq which is U\, satisfies 
the generalized quantum stochastic differential equation. In the derivation of this 
equation an important role plays the quantum dipole noise. 



7. Higher Order Corrections to the Stochastic Limit of the Vacuum 
Expectations 

Let us write some results for vacuum expectation of the evolution operator for the 
models described in previous section. 

In pQ Hamiltonians with polynomial self-interaction have been considered. In 
our case if the system operators D and D + are complex numbers then the expression 
for the vacuum expectation (ABC-formula) has the form 

< U x {t/X 2 ) >=exp( J 4t + A 2 B + A 2 C(t/A 2 )). 

For the linear model, when D and D + are numbers, one has 

A = -jqD+D, B = jiD+D 



Quantum Multipole Noise 



and C(t) -y as t -> oo. Therefore X 2 C{t/X 2 ) = o(A 2 ) and 

< U x (t/X 2 ) >= e-T" tD+D (l + X 2 11 D+D + o(A 2 )). (7.39) 

Now let us consider the case then D and D + are arbitrary, not necessarily 
commutative, system operators. Such models were described in the previous section. 
The expressions for vacuum expectations of the evolution operator for such models 
were obtained in [22] by direct calculations. For the model with rotating wave 
approximation one has 

< U x {t/X 2 ) >= e-i° tD+D [1 + \ 2 jiD + D(l - j tD+D)] 

OO r k 

_ A 2 7l j2 ^ J y.( d+d ) p ~ 1d+2d2 ( d+d ^ p + °( a2 )- ( 7 - 4 °) 

fc=l ' p=l 

In l|7.40|l and (|7.41|) < • > means the averaging only over the boson vacuum so that 
the result is an operator acting in the system Hilbert space . Note that if D and 
D + arc numbers then the series in H7.40J) can be summarized and we obtain fT35jl . 
For the spin-boson model one has 

< U(t/\ 2 ) >„ QC = e lAlt [l - A 2 (Si + Cit)] DD+ 

+ e lA2t [1 - A 2 (B 2 + C 2 t)] D+D + o(A 2 ) (7.41) 
with constants (I = 1,2) 

A, = , dk ^ , B, = [ dk ^ 



wi(fc)-iO' J (wi(Jfc) - iO) 2 ' 

1 



(ui(ki)-i0){u(ki) + u(h2)) 



Lui(ki) — iO cJi(fe) — iO 



Here u>i(k) = u(k) — A, u>2(k) = ui(k) + A > 0. The matrices D and D + in this 
case are 2 x 2-matrices and will be defined in section 9. 

8. The Normally Ordered Form of the Equation for the Evolution Op- 
erator 

Let us write the formal equation for the evolution operator l|6.38|) with the inter- 
action Hamiltonian (|6.31l) using the expansion l|6.32|) of the collective operators. 
Keeping terms of O(X) one has 

= (D(b+(t) + Xc+(t)) - D+(b(t) + Xc(t)))U(X, *) + ... (8.42) 

ot 

Here operators b(t) and c(t) are the white and dipole noise with commutation rela- 
tions JEHU). 

Let us expand the time rescalled evolution operator as the series at the coupling 
constant 

[/(A, t) = U x {t/X 2 ) = U Q (t) + XUi(t) + ... 
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Then comparing terms of the same order of A in the LHS and RHS of equality 
one obtains a system of stochastic differential equations 

jU {t) = (Db + (t) - D+b(t))U (t) (8.43) 

jU x {t) = (Dc+(t) - D+c(t))U (t) + (Db+(t) - D+b^U^t). (8.44) 

The normally ordered form of equations of such kind is the form in which creation 
operators are on the left side of the evolution operator and annihilation operators 
are on the right one. Such form of the equation is convenient for calculation of 
matrix elements of the evolution operator on coherent (exponential) vectors. The 
Eq. Ij8.43jl and Eq. 18.44jl are equivalent to the following 

j t U Q {t) = Db+{t)U Q {t) - D+U Q (t)b(t) - D+[b(t), U (t)} 

j t U x {t) = Dc+(t)U Q (t) - D+U (t)c(t) + Db+^Uxit) - D + U x {t)b{t) 

-D+[c(t),U (t)]-D+[b(t),Ui(t)]. 

Therefore in order to bring the Eq. (|8.43|) and Eq. 18.44f) to the normally or- 
dered form one needs to compute the commutators [b(t),U~o(t)], [c(t),Uo(t)] and 

[&(*). EM*)]- 

The first commutator had been computed in [2] and equal to 

[b(t),U (t)]=j DUo(t) 

The second equal to zero because white noise commutes with the dipole noise and 
Uo(t) "consists" only of the white noise. The last commutator can be computed 
using causal commutation relations JEHU, time consecutive principle and integral 
form of the equation 18.44f> 



b(t),U 1 (t)]=D / dt 1 [b{t),b + (t 1 )]U 1 (t 1 )=^ DU 1 (t). 



Using these commutators one immediately obtains 

jU (t) = Db+(t)U Q (t) - D+U (t)b(t) - l0 D + DU (t) (8.45) 

^-U x {t) = Dc + (t)U a (t) - D+U (t)c(t) 
dt 

+ Db + (t)Ui{t) - D + Ui{t)b{t) - la D+DUi{t). (8.46) 

We can rewrite Eq. I|8.45|l and Eq. I|8.46|l in the forms l|1.8|l and H1.9[) correspond- 
ingly. The equation in the form (|8.46|) we call the generalized quantum stochastic 
differential equation. It is an important task to develop a theory of equations of such 
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form and then to prove the existence and uniqueness of the solution. An interesting 
problem is to study the unitarity of the solution. 

9. Spin-Boson Hamiltonian 

Let us consider the spin-boson model with Hamiltonian i|(j.37fl . The free evolution 
of a z is 

a z (t) = e itHs o z e- ltH ° = e itA D + e - itA D+ 

with 

-H-\ -i 

It is easy to verify that D 2 = D +2 = and (DD + ) 2 = DD + . For this interaction 
(see HJj) 

V(t) = D®W + {t)+D+ ® W{t). 
Here W{t) = A+(t) + A 2 (t) with 

Ai(t) = J dkg(k)e- ltuJl{k) a{k), I = 1,2 

and ui(k) = oj(k) - A, ui 2 (k) = w(k) + A > 0. 

We will use the following expansion for the vacuum expectation value of 
U(\,t) = U x (t/X 2 ): 

< U(X,t) >=< U (t) > +\ 2 < U 2 (t) > +o(A 2 ). 

Theorem 9.3. The first correction f(t) =< U-zit) > satisfies the equation 



dm 

dt 



{iA l DD + + iA 2 D+D)f(t) - DD+de lAlt - D + DC 2 e 



JA 2 t 



(9.47) 



with constants (I — 1,2) 

\g(k)\ 2 



A, 



dk 



u>i{k) — iO ' 



B, 



dk 



\g(k)\ 2 



Q = iAiBi-iZi, 



dk\dk 2 



\g{ki)\ 2 \g(k 2 )f 



{u) l {k 1 )-iQ){L0{k 1 )+^{k 2 )) 



1 



+ 



1 



^i(ki) — iO cui(k 2 ) — iO 



Proof. One has 



a<u ™> =-{< m #mx.t)> 



-i < 
1 



jA 1 (t/\ 2 )(g>D + jA 2 {t/\ 2 )®D+ ) .n\J) 



> 



dt' 



Ai(t/A 2 ),A+(t'/A 2 ) ®DD+ < U(X,t') > 



Quantum Multipole Noise 



+ < A 2 {t' /X^A^t/X 2 ) ® DD + U(X,t') > 



A 2 {t/X 2 ),A+(t'/X 2 ) (g) D + D < U(X, t') > 



+ < A 1 (t'/X 2 )A 2 (t/X 2 ) <g> D+DU(X,t r ) > 



(9.48) 



Let us consider the second term in the first curly bracket. Using the expression 

t 2 t ti 

U{X,t) = 1- i fvih/X^dh + j dh f «ft 2 7(ii/A 2 )V r (t 2 /A 2 )^(A,t 2 ) 



o 







one gets 



lim -i /" < J 4l(</A 2 )A2(^/A 2 )(8)Z?L' + C/(A,^ , ) > 



t ti ta 

= lim^6 / *i / d<2 J dh<A 1 {t/X 2 )A 2 {t 1 /X 2 )W+{t 2 /X 2 )W{h/X 2 ) 



® (£>£>+ ) 2 E/(A,i 3 ) > • 
We suppose that this limit is equal 

t ti ta 

limi /"tfti y dfe y dt 3 <A 1 (t/X 2 )A 2 (t 1 /X 2 )A+(t 2 /X 2 )A+(t 3 /X 2 )> 
ooo 



®{DD+Y < U (t a ) > 

t/X 2 t, t 2 

( J DD+) 2 lim hm/ dfc x y dfc 2 |g(fci)| 2 | 5 (/c 2 )| 2 y dhj dt 2 J dt^e^ " 



x Li^(k 1 )t 3 +tu J2 (k 2 )t 2 +e tu J -(k 2 )t 3 )+tu 2 (k 1 )t 2 } < > = iDD + Zl < f/ (t) > . 



Here we denote (k) := Wi(fc) — ie. 

Using Theorem l4.2l one gets for the first term in the first curly bracket in 



A 2 



Ai (t/X 2 ), Af(t'/X 2 ) (g>DD+ < U(X,t') > dt' 



iA l DD+ < U (t) > +iX 2 A 1 DD + f(t) - X 2 DD + B 1 



d < U (t) > 
dt 



o(A 2 
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Making similar computation with the second curly bracket in l|9.48|l one gets 
d < U f > + X^M = { lAl DD+ + i^D) < U (t) > 

+iA 2 (AiL>Z) + + A 2 D+D)f(i) - \ 2 DD+de lAlt - A 2 D+ DC 2 e lMt ■ 

Let us compare the terms before equal degrees of A in both sides of the last 
relation. For zero order term one has 

d< U °® > = {iA x DD+ + iA 2 D+D) < U (t) >^> 
dt 

< U (t) >= e lAlt DD+ + e iA2t D+D. 
Now we substitute this solution for < Uo(t) > in the terms before A 2 and obtain 



the equation 

dm 

dt 



(iA x DD + + iA 2 D+D)f(t) - DD+C ie lAlt - D + DC 2 e iM 



This finishes the proof of the theorem. □ 
It is easy to verify that a solution of equation l|9.47ll is 

f(t) = [Bx - C 1 t]e iAlt DD+ + [B 2 - C 2 t]e iA ' 2t D + D . 

Here constants Bi, B 2 are determined by initial conditions. In order to obtain (|7.41|) 
one should set B\ = — B\, B 2 = — B 2 . 

10. Conclusions 

In the present paper we have considered the operator valued distribution with com- 
mutation relations proportional to the derivative of 5-function, the so-called dipole 
noise operators. We have constructed a representation of these commutation re- 
lations in a Fock space with indefinite metric, i.e. in a pseudo-Hilbert space. An 
interesting problem is to construct a classical analog of this quantum dipole and, 
more generally, multipole noise. 

The multipole noise is interesting not only as a mathematical object but it also 
plays an important role for applications in physics. In particular, we apply the dipole 
noise operators to study some models of quantum optics. For the evolution operator 
of these models we obtain stochastic equations which we call generalized quantum 
stochastic differential equations. An important open problem is to develop a theory 
of these generalized quantum stochastic equations, including proving of existence 
and uniqueness of the solution. An interesting question is to study an analogue of 
the unitarity of the solution. 

The multipole noise describes higher order corrections to the stochastic limit of 
quantum theory. This is an approximate non-perturbation method which one can 
apply to study the dynamics of quantum system interacting with a reservoir in the 
case when coupling constant or density of particles of reservoir is a small parameter 
and for large time. Therefore it is an important task to apply this approach to other 
models of quantum theory. 
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